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The main goal of this paper is deriving Density of states g{t) (degeneracy function) 
per volume for an equation of state (EOS) p = —p (we called it dark energy(DE)). We 
have concluded that thermodynamic quantities such as pressure and energy density 
, are simple functions of temperature, fugacity, curvature and mass of Bosons. Our 

work has been expressed the origin of some claims about the negativity of the entropy 
for the scalar fields models of DE. 
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I: INTRODUCTION 



In recent years some works have been done on the thermal properties of the negative 
pressure matters, especially on the dark energy. The main assumption of all these models 
is that the DE is a thermal ensemble at a certain temperature with a specific amount of 
entropy [l]]. It would be usable to take temperature as an intrinsic characteristics of DE, 
namely a barotropic matter with an EOS as p = f(p). This fluid must not violate energy 
conditions. Also we can take both dark matter (DM) and DE as different types of the same 
matter field with a phase transition in cosmo logical distances There is not any reason- 
able statistical description in negative pressure matters .In some of recent papers, thermal 
properties of dark energy have been discussed in terms of the assumption that the dark 
energy substance is a thermal ensemble at a certain temperature with an associated thermo 
dynamical entropy [1]. It is usually assumed that this temperature is an intrinsic property 
of DE rather than the temperature of the heat bath fixed by surrounding matter. 
Checking the second law of thermodynamics in the context of the black hole physics is 
another important problem. Also the generalized second law (GSL) must be satisfied by 
any physical description of the DE even in the new scenarios of the modified gravity in 
the non relativistic diffeomorphism broken models such as Horava proposal [29]. If we de- 
note that S to t is the total entropy of a system,then GSL guaranties that S to t > Ofor all 
times. Recently Mubasher Jamil , et.al have shown that [26] in an interactive model for DE , 
DM and radiation " the generalized second law is always and generally valid, independently 
of the specific interaction form of the fluid's equation- of- state parameters and of the back- 
ground geometry." .The history of GSL backs to the Unruh and Wald classical work[28]and 
Bekenstein works on the black holes [31]. There exists a recondite relation between the holo- 
graphic dark energy model and generalized second law of thermodynamics as in [27] and 
also in quintom dominated universe [30]. Holographic scenario for DE is so popular and was 
investigated by authors both in the context of the usual FRWL cosmology and also in the 
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modified ones and in a braneworld picture of universe[25,32,33,36].There is a delicate re- 
lation between Bulk brane interaction , holographic dark energy [34], Gauss Bonnet dark 
energy models[35, 38], Stringy inspired tachyon model [37], phantom like regime [39], Observa- 
tional constraints on holographic dark energy [40], Holographic Chaplygin gas model[41], 
Holographic Chaplygin DGP Cosmologies [ 42] a nd even Holographic Modified Gravity [43]. 

Some years ago Kulikov and Pronin jsf, THJ constructed a simple formulation of a local 
quantum statistics of Bosonic field in a curved background. Indeed they explained a simple 
expression for a grand canonical thermodynamic potential density Q describing locally all 
thermo chemical properties of gases with a curvature dependent parts [5]. Later The high- 
temperature expansion of the grand thermodynamic potential of non conformably invariant 
spin-0 and spin-1/2 gases in an arbitrary static spacetime with their boundary is calculated 
through the method conforming from Kirsten [16]. 

In this report, we have worked on the statistics of a negative matter with EOS like DE. At 
first we have obtained the degeneracy distribution function for massive Bosons in this grand 
ensemble and it has been shown that the density of states for such an exotic matter is a 
constant function of energy . We have summed over all energy states once more and obtained 
a closed form for the grand canonical thermodynamic potential density and then we have 
inferred the expressions for the density of the entropy and other important functions. We 
have also shown that these thermodynamic equations are generally related with the fugacity , 
the Ricci scalar and consequently the temperature and the chemical potential. The chemical 
potential has been associated with a conserved Boson number. We have stated that from 
a statistical point of view DE violate energy conditions companionship to the claims stated 
in the recent papers^. To shunning from this undesirable result we must take both volume 
and temperature as two dependent variables as Gong and collaborations have shown in [6]. 



II: LAGRANGIAN OF THE MODEL AND SOME REMARKS ABOUT THE 

BACKGROUND METRIC 

We begin from Lagrangian of a massive scalar field coupled to a gravity which we assume 
that is non negative. Consider the action js|] 

S = J v^ja - \(m 2 + iR)^ 2 ) (1) 

where R is the scalar curvature and we assume that the metric of the space time g^ u is static 
or slowly time variable such that we can take the Ricci scalar, time independent. Spherical 
symmetry is not essential in this formalism, because we use the general form for the metric, 
and only we impose a restriction on the time variation of the metric. The only essential 
assumption is on the surface tension energy portion in EOS. In the case of bulk matter we 
ignore the surface effects and all of the thermodynamics quantities are expressed per unit 
volume. Considering the equilibrium assumption , a static metric is essential [14]. In this 
case the background metric g^ v generated by the mass distribution is assumed to be a slowly 
varying function of the time on the inverse temperature scale[14]. More precisely, we assume 

at << p ' 1 } 

The fluctuation of the metric is due to the energy momentum tensor fluctuations. If we ignore 
the energy fluctuations, (equilibrium state) , we can assume that the metric is static or slow 
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varying with respected to the time. We neglect the influence of the matter and radiation and 
also assume that their interactions with DE are small and serve only to provide a heat bath 
at the temperature T. The boundary is a timelike tube which is periodically identified in the 
imaginary time direction with the period (3. The necessary path integral for the construction 
of the partition function is taken over asymptotically vanishing fields which are periodic in 
the imaginary time t with the period (3. Thus, the functional integration has been assumed 
with the periodicity in the imaginary time and the asymptotic flatness of the metric fields. 
As we know that in any static or slowly time varying spacetimes the Ricci scalar has no 
dynamics (time dependency) or it's variations with respect to the time is negligible, but only 
in these models we can write a suitable explicit expression for Green function. On the other 
hand the static spacetime implies that the ensemble of Bosons remains at a definite constant 
temperature in an equilibrium or a quasi equilibrium state of the system .If we are far away 
from a static space times, all equations has only a perturbative solution and it's applications 
of an expanding universe is doubtful. There is only some few published works in this topics 
and finally leads to a stochastic statistical field theory [20]. 1 The topology of this manifold 
is considered to be a S* 1 x M 3 and also we can take M 3 manifold without any boundary [17] 



III: THERMODYNAMIC POTENTIAL 

In this section we will analyze the scalar field model with a conserved charge. 
The Lagrangian of the model is: 

S m = J d A x^gJx)^{x){-U x + m 2 + ZR)Q(x) (3) 

Where $ = (0i,02) is a doublet of the real fields. The action written in the terms of real 
fields will be: 

S m = ~{\) J d 4 x^/gJx)<p a (x)(-D x + m 2 + ZR)4> a {x) (4) 

The total action of the system "matter + gravity" is: 

Stot — S g + S m (5) 
Now we can write the effective action at finite temperature as: 

L eff (P) = L g - U tf,n,R) (6) 



Stochastic semi classical gravity in the 1990s is a theory naturally evolved from semi classical gravity in 
the 1970s and 1980s. In stochastic semiclassical gravity the main object of interest are the noise kernel, the 
vacuum expectation value of the (operator-valued) stress-energy bi-tensor, and the centerpieces being the 
(semiclassical) Einstein-Langevin equation. It also brings out the open system concepts and the statistical 
and stochastic contents of the theory such as dissipation, fluctuations, noise and decoherencc. Hu[20] has 
been described the applications of the stochastic gravity to the back reaction problems in the cosmology 
and the black-hole physics. Further discussions of the ideas and the ongoing research topics can be found 
in [21,22,23] 
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where L g is 



Lg — L g 



■ roc 

- / dm 2 trG S D(x,x) 

Z Jm? 

The symbol tr(...) is determined as: 



and u;(/3, i?) is the density of grand thermodynamic potential. 

The result (6) may be obtained with the momentum space representation for the Green's 
function of a Boson . 

In the momentum space representation, the expression for L e ff is splited into two parts: 

• poo 

Atf/ = -y / dm 2 trG(x,x)-cj(P,n,R) (8) 

Z J m 2 

The potential u(/3,R) is: 

I fOO 2 OO rs n. ,3, 

= -(-)tr dm^^m-Q^y x / (^3^ + ^ ( 9 ) 

2 



= (^)E7i(i2)(-^*rM^ + 6 2 ) 

where u n = and the geometrical coefficients "fj(R) can be represented in terms of 
scalar R and coupling constant £ as: 

lo(R) = 1 

For introducing the chemical potential we will change the Matsubara frequencies [24] 
oo n — > w n + /x and then thermodynamic potential will be u)(/3, n, R). 
Since both positive and negative frequencies are summed, we will get: 

tr ln(w 2 + e 2 ) ->■ ir ln[(w n + /i) 2 + e 2 ] 
= tr{\n[ul + (e - /x) 2 ] + ln[u, 2 + (e + /x) 2 ]} (10) 

After doing required mode's summation in (9) we will have: 

u(f3, n, R) = fi, R) + u + ((3, R) (11) 

Where: 



(A K R) = (h £ 7i(i?)(- Ay ln(l - exp[-/3(e - /,)]) (12) 



And: 

2 



u+(P,KR) = (^)E7,( J R)(-^ln(l-^exp[-/3(e -/,)]) (13) 

p i=o 

So the density of grand thermodynamic potential is the series: 

2 

u>{p,n,R) = Y,n(.R)h{P™,z) ( 14 ) 

i=o 

Where: 

b (pm,z) = (^)ln(l - zexp(-/9e)); 

b j (Pm,z) = (--£-yb (Pm,z) (15) 
and the fugacity is z = exp(/fyt). 

IV: STATISTICS AND THERMODYNAMICS OF AN IDEAL (NON 

INTERACTIVE) BOSE GAS 

The Bose distribution function as the derivative of the grand thermodynamic potential 
is given by: 



d/j, 



ns = _ ™ k KW,~J (16) 



For occupation numbers with momentum k we can obtain: 

1 

(z-ie&k - 1) 

Where the function B(/3,R) is described by the formula: 



"£= TITOi— -, B{P,R) (17) 



Btf, R) = l+ li(R)^-[l ~ (1 - ^exp-^)- 1 ] + ... (18) 

The function B(j3, R) depends on the curvature, the temperature and the energy of the 
Boson[18]. 

Studying the thermodynamic properties of the Bose gases we will start with the equation: 



iP^R) = -(\)j^iAR){-^) j W-^M-^\) (19) 



d 
om 

3=0 



In the non-relativistic limit for the particle energy e = we can derive from (19) the 
equation: 

2 r, 

u(P,vl,R) = ^ 7i (^ 6/2 ( 2! )(-_)iA- 3 (20) 

3=0 
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where A = (2n /mT) 1 / 2 is a wavelength of the particle, and the function g 5 / 2 (z) has the 
following form: 



00 4 f 



00 

#5/2(2) = Vy^ J= I x 2 \n(l - z exp(-x 2 ))dx 



The average number of the particles in a certain momentum state k is obtained as the 
derivative: 

< n t >= —q-u(P,1*,R) = ^7,(^)# 5 /2(2)(-^) j (^ 1 exp(^e) - 1) 

The density of the particles is: 

n = A- 3 [l - 7l (i?)(3/4m 2 ) - 72 ( J R)(3/16m 4 )]( ? 3/2(2) + n (21) 
Where the new function g 3 / 2 (z)is: 

9s/2(z) = z^-g 5/2 (z) 

and 

z 



n = 



l-z 



is the average number of the particles with zero momentum. The functions #3/2(2) and 
#5/2 ( z ) are special cases of a more general class of functions: 

00 

= X>V* B (22) 
fc=i 

In a more simple form the equation (21) may be written as: 

(n-n )\ 3 = g 3/2 {z,R) (23) 

Where: 

g 3/2 (z,R) = [l-a^ + ...]g 3/2 (z) (24) 

is a function which depends on the curvature, and a is a numerical parameter. 

The equation (23) connects four values: fugacity, temperature, density of the particles and 

curvature. 

V: CALCULATING THE DENSITY OF STATES FOR EOS p = -p 

The density of grand thermodynamic potential is : 
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2 

u(J3,p,R) = y - i^ 7i (i2)(-^)i^ln(l - zexp(-Pe)), (25) 
where = 1/T. 

We can write the sum over the energy states as an integration over energy values by 
considering it's occupation number function (Degeneracy of the energy levels) g(e) where: 

V 

We know that for photons, phonons and non-relativistic particles with spin s this function 
is(in Geometrical units where: h = c = 1) [ll| : 

Ve 2 
(2tt) 
9Ne 2 

g(e) = — —{phonons) 
u D 

, , (2s + l)\/m 3 /V/% Ar , . . . „ 
g{e) = ^_ {J\ onrelativisticBosegas) 

We now review some relations between the Grand Thermodynamic Potential of the system 
Q((3,z,V), the Helmholtz Free energy (F), the entropy (S)and the pressure (P) [ill ]: 

F = U - TS 

U = Vl + fjiN + TS 

pv = -n 
b ~ of 

Q = p(Ln(Z(fi,T, TV)) 
dLn(Z{fi,T,N)) 
1 d(3 J ' 

where in these relations z is called fugacity : 

z = e p » 

Since now we named the density of the Grand potential of the system as u then we can 
write next equations for the density ( per volume) of them as: 

f = p-Ts (26) 
p = u + fin + Ts (27) 
P = -co (28) 

The main goal of this paper is deriving degeneracy function per volume for an EOS p = —p, 
and then if it is the well behavior function, we can fit it with an interaction inside a Boson 
gas, then the generation function of both systems will become the same. Our result for 
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the dispersion relation has been shown that there is the k 3 dependence of the energy. This 
dependence can be the main effective term of the interaction in a definite density and 
temperature, and also like the logarithmic dependence of the dispersion relation strength 
quark matter (SQM); the k 3 term can be dominant in the some regions[19]. 

Remembering previous relations together with EOS p = — p we lead to the next first 
order partial differential equation for u. 

du n du: 

= p w (29) 

now we note here that the general solution for above equation in which 

u( l l,P) = Y,CaW (3°) 
a 

we have a simple boundary condition which is concluded from thermodynamics: the entropy 
of a system must be vanish at absolute zero temperature (3. that is we must have: 

lim (A = 

/3^oo dp 

Thus we must limit ourselves only to a < —1. But the coefficients remain undetermined. 
Thus we take another endeavor to solve this equation. We substitute ( 25) in (29) and take 
derivatives and equalize from both sides all terms which have the same order of m 2 (i.e.; the 
coefficients of (— , j = 0, 1, 2 We have: 



Jo 



e #( e ) 1 



i r°° 

-de-- ln(l - ze-^)g{e)de (31) 
1 P Jo 



We assume that g{e) = an d also we suppose that <r(e) is a well behavior function of 
energy e at the ground state energy level e = 0. After integrating part by part and expanding 
both sides in terms of z (0 < z < 1) and set equal all terms of same z power we have: 



oo 



e-^[a(e) - Ade = 4^, n > 1 (32) 
de pn 

We can define Laplace transformation of the function <r(e) as a function of parameter s = (3n 
by: 

POD 

S(s) = / e~ se a(e)de (33) 
Jo 

We arrive at the next differential equation for the transformed function: 

Z(s) + %-(sX(s) - v(0)) = — (34) 

OS s 



Which has the solution: 



= — - + — (35) 
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In this function C\ is an unknown constant that can be determined using initial conditions. 
But we take it as a constant and find it's value some later. Because we solve the partial 
differential equation (29), then the result can be a function of other quantities like the mass 
of particles. The mass is not an independent thermodynamic variable. If we assume the mass 
as a thermodynamic variable, then the derivative of mass with respected to the density, the 
temperature or the chemical potential must be occurred in the fundamental thermodynamic 
relations [19]. 

By taking inverse Laplace transformation from (35)we obtain: 

(7(e) = (7(0) + Cie (36) 

We mention here that for a Boson's system the chemical potential is always negative namely 
— oo < /i < 0. Since g(e) = we finally obtain: 

9(e) = C, (37) 
Substituting (37) in definition of g(e) 2 we lead to the following dispersion relation for DE: 

k 3 = 6n 2 de - 6tt 2 C 2 (38) 

We know that this is not a dispersion relation for a non interactive matter field. 

But there is no simple system which has such distribution. If we assume that the energy 
levels of such Bosonic system must be positive and the momentum of any particle also be 
positive, then the energy levels of this system must satisfy the following inequality: 

6 " % (39) 

Without any dismay from loss of generality we can take C 2 = .With this choice we accept 
that our ground state energy level is located at e = 0. Another pickup is tagging this energy 
level a non zero value. This situation could be solved by a shift of any energy value up to 
this value. 

Substituting (37) in (25) and making use from: 

°° z n 
1 

n=l 

finally we obtain next function for the density of grand thermodynamic potential: 

1 2 r) 

u(vl,R,I3) = -^Polylog(2,z)J2^(R)(-Q-^y Cl ^ 

p j=o m 

Now we calculate the density of the average of the energy We can label it as pa which is: 



E r 9(e)e , 
PA= V = J ^ 



2 Note that V is disappeared since we calculate density. 
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After substituting (37) in it and doing the integration we obtain: 

p A = j±Polylog(2,z) (41) 
The mean value of the number of the particles in the unit volume of the system is: 

Since < z < 1 the value of n is not negative. Using s = y = — |^ we have: 

2 

s = --Poly\og{2,z)Y j a 3 {R){-^YC l (43) 

' 3=0 

Because the mass dependence of g(e) is undetermined , so we can not conclude strongly 
about the sign of entropy. It can be negative or positive [14]. 

Negative entropy is problematic if we accept that the entropy is in an the association 
with the measure of the number of microstates in statistical mechanics. The intuition of 
the statistical mechanics requires that the entropy of all physical components have to be 
positive. Besides if we consider the universe as a thermodynamic system, the total entropy 
of the universe including DE and DM should satisfy the second law of thermodynamics. 
The GSL for phantom and non-phantom DE has been explored in [12]. It was found that 
the GSL can be protected in the universe with DE. The GSL of the universe with DE has 
been investigated in [13,25,26,27] as well. 3 ,and also the equation of state of the DE is 
uniquely determined and the phantom entropy is negative [13]. For the phantom,when we 
running to the negative entropy problem, the GSL is violated. 4 Considering the apparent 
horizon as the physical boundary of the universe, it was found that both the temperature 
and the entropy can be positive for the DE, including phantom. In the general framework 
that we discussed, where the cosmological metric has a slow varying time dependency, we 
can assume that the time contiguity of the volume and the temperature is venial and for this 
reason we have obtained a negative entropy. We must note here that if in an ensemble of 
particles the volume and temperature were non static, the statistical behavior of this system 
is not in the context of the equilibrium thermodynamics. 5 



VI: CONCLUSIONS 

We have derived a grand canonical partition function for a DE's like EOS. The thermo- 
dynamic equations are generally expressed in terms of three variables: the temperature T , 



3 We may think the DE temperature is equal or proportional to the horizon temperature Tjj- 

4 It is more realistic to consider that the physical volume and the temperature of the universe are related, 
since in the general situation they both depend on the scale factor a(t).It was found that the apparent 
horizon is a good boundary for keeping thermodynamics' laws [13]. 

5 As in a FRW standard model of cosmology, since our background metric is not static and thereupon 
the apparent horizon is time variable, the construction of a local quantum statistical and counting all 
microstates must affect the negative entropy and we can deduce a statistical proof for non negativity and 
non violation of GSL like Gong and collaborations proof [6] . 
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chemical potential /i associated with a conserved particle number N, and scalar curvature R. 
We also reviewed some known results for an ideal Bose gas in curved spacetime. Comparing 
our results with an ideal gas, we have concluded that there is a significance generic difference 
between the statistical distribution of negative pressure matters and ideal gases. Specifically 
, the dispersion relation for a dark energy like EOS is non quadratic (in contradiction with 
a non interactive ideal Bose gas) and comparing this dispersion function with the statistical 
model of SQM we observe that also in SQM systems this dispersion relation is valid in some 
physically accessible region of the system [19]. 

Also we have derived the mean value of the energy of the system in terms of the curvature 
. We have shown that if a matter field (Boson) obeys an EOS like DE, and in cosmological 
language violates the well known energy conditions in General relativity , the degeneracy 
function of this gas must be an implicit constant function of the energy. This function 
implicitly is dependent on the mass of the particles which is not a thermodynamic quantity 
and remains as an undetermined function in this formula. 

The derived thermodynamic equations has been shown that the entropy may be positive 
or negative. As it was stated in [6], this negativity of the entropy is hidden behind this 
significance assumption that the volume and the temperature are independent. But we 
know that in an expanding FRW universe both of them are functions of time. Hence if we 
assume that our metric is quasi static or has a negligible dynamics, any phantom field of 
DE with an intelligible statistical testimony has a negative entropy. If we want to generalize 
our argument to a realistic cosmological case, we must work in a non equilibrium statistical 
regime. 
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found the corresponding entropy and also our work exposure the statistical pedigree 
of some claims about scalar fields which presenting negative pressure. 
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3 1 - I: INTRODUCTION 

CT 

In recent years some works have been done on the thermal properties of negative pres- 
sure matters specially on the dark energy . The main assumption of all these models is 
that the DE is a thermalized ensemble at a certain temperature with a specific amount of 
entropy [3]. It would be usable to take temperature as an intrinsic characteristics of DE, 
namely a barotropic matter with an EOS as p = f(p). This fluid must not violate energy 
conditions. Also we can take both dark matter (DM) and dark energy as different types 
fSJ ■ of the same matter field with a phase transition in cosmological distances 0]. There is not 
^ . any reasonable statistical description in negative pressure matters .In some of recent papers, 
thermal properties of dark energy have been discussed in terms of the assumption that the 
dark energy substance is a thermalized ensemble at a certain temperature with an associ- 
ated thermo dynamical entropy [1]. It is usually assumed that this temperature is an intrinsic 
property of DE rather than the temperature of the heat bath fixed by surrounding matter. 
Checking the second law of thermodynamics in the context of the black hole physics is an- 
other important problem.Also the generalized second law (GSL) must be satisfied by any 
physical description of the DE even in the new scenarios of the modified gravity in the non 
relativistic , diffeomorphism broken modeles such as Horava proposal [29]. If we denote that 
S to t is the total entropy of a system,then GSL guaranties that S to t > Ofor all times. Recently 
Mubasher Jamil , et.al showed that [26] in an interactive model for DE , DM and radia- 
tion " the generalized second law is always and generally valid, independently of the specific 
interaction form of the fluid's equation- of- state parameters and of the background geome- 
try." .The ancestry of GSL backs to the Unruh and Wald classical work[28]and Bekenstein 
works on the black holes [31]. There exists a recondite relation between the holographic dark 
energy model and generalized second law of thermodynamics as was pointed out in [27] and 
also in in quintom dominated universe [30]. Holographic scenario for DE is so popular and 
was investigated by authors both in the context of the usual FRWL cosmology and also in 
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the modified ones and in a braneworld picture of universe[25,32,33,36].There is a delicate 
relation between Bulk brane interaction , holographic dark energy [34], Gauss Bonnet dark 
energy models[35, 38], Stringy inspired tachyon model [37], phantom like regime [39], Observa- 
tional constraints on holographic dark energy [40], Holographic Chaplygin gas model[41], 
Holographic Chaplygin DGP Cosmologies [ 42] a nd even Holographic Modified Gravity [43]. 

Some years ago Kulikov and Pronin [3|, |l5j constructed a simple formulation of a local 
quantum statistics of Bosonic field in a curved background. Indeed they explained a simple 
expression for a grand canonical thermodynamic potential density |4J which locally describes 
all thermo chemical properties of gases with a curvature dependent parts [EJ. Later The high- 
temperature expansion of the grand thermodynamic potential of non conformably invariant 
spin-0 and spin-1/2 gases in an arbitrary static spacetime with their boundary is calculated 
by method conforming from Kirsten [16]. 

Through this report, we have worked on the statistics of a negative matter with EOS like 
DE. At first we have obtained the degeneracy distribution function for massive Bosons in 
this grand ensemble and it has been shown that the density of states for such exotic matter 
is a constant function of energy . We have summed over all energy states once more and 
obtained a closed form for grand canonical thermodynamic potential density and then we 
have inferred the expressions for the density of entropy and other important functions. We 
have also show that these thermodynamic equations are generally expressed in terms of 
fugacity , Ricci scalar and consequently on the temperature and chemical potential. The 
chemical potential has been associated with a conserved Boson number . We have been 
stated that from a statistical point of view DE violate energy conditions companionship 
to the claims stated in the recent papers jH]. To shunning from this undesirable result we 
must take volume and temperature as two dependence variables as Gong and collaborations 
showed in [6]. 



II: LAGRANGIAN OF THE MODEL AND SOME REMARKS ABOUT THE 

BACKGROUND METRIC 

We begin from Lagrangian of a massive scalar field coupled to gravity which we assume 
that is non negative. Consider the action js|] 

8 = j >/=£d*s(i^ - i(m a + iR)f) (1) 

where R is the scalar curvature and we assume that the metric of space time g^ u is static 
or slowly time variable such that we can take the Ricci scalar, time independent. Spherical 
symmetry is not essential in this formalism, because we use the general form for the metric, 
and only we impose a restriction on the time variation of the metric. The only essential 
assumption is on the surface tension energy portion in EOS. In the case of bulk matter we 
ignore the surface effects and all of the thermodynamics quantities are expressed per unit 
volume. Considering the equilibrium assumption , a static metric is essential [14]. In this 
case the background metric g^ u generated by the mass distribution is assumed to be a slowly 
varying function of the time on the inverse temperature scale[14]. More precisely, we assume 



at p 
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The fluctuation of the metric is due to the energy momentum tensor fluctuations. If we ignore 
the energy fluctuations, (equilibrium state) , we can assume that the metric is static or slow 
varying with respected to the time. We neglect the influence of the matter and radiation 
and also assume that their interaction with DE are small and serve only to provide a heat 
bath at the temperature T. The boundary is a timelike tube which is periodically identified 
in the imaginary time direction with the period /3. The necessary path integral for the 
construction of the partition function is taken over asymptotically vanishing fields which 
are periodic in imaginary time t with period (3. Thus, the functional integration has been 
assumed with the periodicity in imaginary time and the asymptotic flatness of the metric 
fields. As we know that in any static or slowly time varying spacetimes the Ricci scalar has 
no dynamics (time dependency) or it's variations with respect to the time is negligible, but 
only in these models we can write a suitable explicit expression for Green function. On the 
other hand the static spacetime implies that the ensemble of Bosons remains at a definite 
constant temperature in an equilibrium or a quasi equilibrium state of the system .If we are 
far away from a static space times, all equations has only a perturbative solution and it's 
applications of an expanding universe is doubtful. There is only some few published works 
in this topics and finally leads to a stochastic statistical field theory [20]. 1 The topology of 
this manifold is considered to be a 5 1 x M 3 and also we can take M 3 manifold without any 
boundary [17] . 

Ill: THERMODYNAMIC POTENTIAL 

In this section we will analyze the scalar field model with a conserved charge. 
The Lagrangian of the model is: 

S m = ~{\) f d 4 xy/gJx)^*(x)(-D x + m 2 + ZR)Q(x) (3) 

Where $ = (0i,02) is a doublet of the real fields. The action written in the terms of real 
fields will be: 

Srn = ~{\) f d 4 X^/^(x)r(x)(-D x + W 2 + 6R)0aOr) (4) 

The total action of the system "matter + gravity" is: 

Stot = S g + S m (5) 
Now we can write the effective action at finite temperature as: 

Leff(/3) — L g — u}(/3, ji, R) (6) 



Stochastic semi classical gravity in the 1990s is a theory naturally evolved from semi classical gravity in 
the 1970s and 1980s. In stochastic semiclassical gravity the main object of interest are the noise kernel, the 
vacuum expectation value of the (operator-valued) stress-energy bi-tensor, and the centerpieces being the 
(semiclassical) Einstein-Langevin equation. It also brings out the open system concepts and the statistical 
and stochastic contents of the theory such as dissipation, fluctuations, noise and decoherencc. Hu[20] has 
been described the applications of the stochastic gravity to the back reaction problems in the cosmology 
and the black-hole physics. Further discussions of the ideas and the ongoing research topics can be found 
in [21,22,23] 
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where L g is 



Lg — L g 



■ roc 

- / dm 2 trG S D(x,x) 

Z Jm? 

The symbol tr(...) is determined as: 



and u;(/3, i?) is the density of grand thermodynamic potential. 

The result (6) may be obtained with the momentum space representation for the Green's 
function of a Boson . 

In the momentum space representation, the expression for L e ff is splited into two parts: 

• poo 

Atf/ = -y / dm 2 trG(x,x)-cj(P,n,R) (8) 

Z J m 2 

The potential u(/3,R) is: 

I fOO 2 OO rs n. ,3, 

= -(-)tr dm^^m-Q^y x / (^3^ + ^ ( 9 ) 

2 



= (^)E7i(i2)(-^*rM^ + 6 2 ) 

where u n = and the geometrical coefficients "fj(R) can be represented in terms of 
scalar R and coupling constant £ as: 

lo(R) = 1 

For introducing the chemical potential we will change the Matsubara frequencies [24] 
oo n — > w n + /x and then thermodynamic potential will be u)(/3, n, R). 
Since both positive and negative frequencies are summed, we will get: 

tr ln(w 2 + e 2 ) ->■ ir ln[(w n + /i) 2 + e 2 ] 
= tr{\n[ul + (e - /x) 2 ] + ln[u, 2 + (e + /x) 2 ]} (10) 

After doing required mode's summation in (9) we get: 

u(f3, n, R) = fi, R) + u + ((3, R) (11) 

Where: 



(A K R) = (h £ 7i(i?)(- Ay ln(l - exp[-/3(e - /,)]) (12) 



And: 

2 



u+(P,KR) = (^)E7,( J R)(-^ln(l-^exp[-/3(e -/,)]) (13) 

p i=o 

So the density of grand thermodynamic potential is the series: 

2 

u>{p,n,R) = Y,n(.R)h{P™,z) ( 14 ) 

i=o 

Where: 

b (pm,z) = (^)ln(l - zexp(-/9e)); 

b j (Pm,z) = (--£-yb (Pm,z) (15) 
and the fugacity is z = exp(/fyt). 

IV: STATISTICS AND THERMODYNAMICS OF AN IDEAL (NON 

INTERACTIVE) BOSE GAS 

The Bose distribution function as the derivative of the grand thermodynamic potential 
is given by: 



d/j, 



ns = _ ™ k KW,~J (16) 



For occupation numbers with momentum k we can obtain: 

1 

(z-ie&k - 1) 

Where the function B(/3,R) is described by the formula: 



"£= TITOi— -, B{P,R) (17) 



Btf, R) = l+ li(R)^-[l ~ (1 - ^exp-^)- 1 ] + ... (18) 

The function B(j3, R) depends on the curvature, the temperature and the energy of the 
Boson[18]. 

Studying the thermodynamic properties of the Bose gases we will start with the equation: 



iP^R) = -(\)j^iAR){-^) j W-^M-^\) (19) 



d 
om 

3=0 



In the non-relativistic limit for the particle energy e = we can derive from (19) the 
equation: 

2 r, 

u(P,vl,R) = ^ 7i (^ 6/2 ( 2! )(-_)iA- 3 (20) 

3=0 
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where A = (2n /mT) 1 / 2 is a wavelength of the particle, and the function g 5 / 2 (z) has the 
following form: 

00 z l 4 f°° 

#5/2(2) = ^2 ]H72 = ~7/^ ^ ^ ln ( 1 ~ 2 exp(-:r 2 ))cfe 

The average number of the particles in a certain momentum state k is obtained as the 
derivative: 

^ j=o m 

The density of the particles is: 

n = A- 3 [l - 7l (i?)(3/4m 2 ) - 72 ( J R)(3/16m 4 )]( ? 3/2(2) + n (21) 
Where the new function g 3 / 2 (z)is: 

9s/2(z) = z^-g 5/2 (z) 

and 

z 



n = 



l-z 



is the average number of the particles with zero momentum. The functions (73/2(2) and 
95/2(2) are special cases of a more general class of functions: 

00 

= X>V* B (22) 
fe=i 

In a more simple form the equation (21) may be written as: 

(n-n )\ 3 = g 3/2 {z,R) (23) 

Where: 

g 3/2 (z,R) = [l-a^ + ...]g 3/2 (z) (24) 

is a function which depends on the curvature, and a is a numerical parameter. 

The equation (23) connects four values: fugacity, temperature, density of the particles and 

curvature. 

V: CALCULATING THE DENSITY OF STATES FOR EOS p = -p 

The density of grand thermodynamic potential is : 



7 



2 

u>(J3,v,R) = y= -E^^^^EM 1 -^-^)), (25) 

P 3=0 

where (5 = 1/T. 

We can write the sum over energy states as an integration over energy values by consid- 
ering it's occupation number function (Degeneracy of energy levels) g(e) where: 

V 

We know that for photons, phonons and non-relativistic particles with spin s this function 
is(in Geometrical units where: h = c = 1) [ll| : 

Ve 2 
(2tt) 
9Ne 2 

g(e) = — —{phonons) 
u D 

, , (2s + l)\/m 3 /V/% Ar , . . . „ 
g{e) = ^_ (iv onrelativisticBosegas) 

We now review some relations between the Grand Thermodynamic Potential of the s yste m 
£l(/3,z,V) and the Helmholtz Free energy of system(F), entropy (S)and pressure (P) [111 ]: 

F = U — TS 

U = Vl + fiN + TS 

pv = -n 
b ~ of 

Q = p(Ln(Z(fi,T, N)) 
dLn(Z{fi,T,N)) 
1 d(3 J ' 

where in these relations z is called fugacity : 

z = e p » 

Since now we named the density of the Grand potential of the system as u then we can 
write next equations for the density ( per volume) of them as: 

f = p-Ts (26) 
p = u + fin + Ts (27) 
P = -co (28) 

The main goal of this paper is deriving degeneracy function per volume for an EOS p = —p, 
and then if it is the well behavior function, we will fit it with an interaction inside a Boson 
gas, then the generation function of both systems will become the same. Our result for 
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dispersion relation shows that there is the k 3 dependence of the energy. This dependence 
can be the main effective term of the interaction in a definite density and temperature and 
like the logarithmic dependence of the dispersion relation strength quark matter (SQM)(in 
the some regions) the k 3 term can be dominant [19]. 

Remembering previous relations together with EOS p = — p we lead to the next first 
order partial differential equation for u. 

"W = "«r (29) 

now we note here that the general solution for above equation is: 

u(ll,P) = Y,CaWY (3°) 

a 

we have a simple boundary condition which is concluded from thermodynamics: the entropy 
of a system must be vanish at absolute zero temperature (3. that is we must have: 

lim (A = 

fS^oo dp 

Thus we must limit ourselves only to a < —1. But the coefficients remain undetermined. 
Thus we take another endeavor to solve this equation. We substitute ( 33) in (37) and take 
derivatives and equalize from both sides all terms which have the same order of m 2 (i.e.; the 
coefficients of (—^Y^j = 0, 1, 2 We have: 

[70=l it = -j[ Ml - ze ^ )gii)di (31) 

We assume that g(e) = an d also we suppose that a(e) is a well behavior function of 
energy e at the ground state energy level e = 0. After integrating part by part and expanding 
both sides in terms of z (0 < z < 1) and set equal all terms of same z power we have: 



e-^W)-Ade= a M n > l (32) 
de pn 



o 



We define Laplace transform of function er(e) as a function of parameter s = (3n by: 

P oo 

S(s) = / e~ S€ a(e)de (33) 
Jo 

We arrive at the next differential equation for transformed function: 

Z(s) + £-(8X(s) - a(0)) = — (34) 

OS s 

Which has solution: 

S( s ) = — - + — (35) 

In this function C\ is an unknown constant that can be determined using initial condi- 
tions. But we take it as a constant and find it's value some later. Because we solve the 
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partial differential equation (29), then the result can be a function of other quantities like 
mass of particles. The mass is not an independent thermodynamic variable. If we assume 
the mass as a thermodynamic variable, then the derivative of mass with respected to the 
density, temperature or chemical potential must be occurred in fundamental thermodynamic 
relations [19]. 

By taking inverse Laplace transformation from (35)we obtain: 



a(e) = <t(0) + de (36) 
We mention here that for a Boson's system the chemical potential is always negative namely 

Me 
de 

9(e) = d (37) 



— oo < /i < 0. Since g(e) = we finally obtain: 



This relation is some eccentric. Substituting (37) in definition of g(e) 2 we lead to the 
following dispersion relation for DE: 

k 3 = §T?C x e - 6n 2 C 2 (38) 

We know that this is not a dispersion relation for a non interactive matter field. 

But there is no simple system which has such distribution. If we assume that the energy 
levels of such Bosonic system must be positive and the momentum of any particle also be 
positive, then the energy levels of this system must satisfy the following inequality: 

6 " % (39) 

Without any dismay from loss of generality we can take C 2 = .With this chooses we accept 
that our ground state energy level is located at e = 0. Another pickup is tagging this energy 
level a non zero value. This situation could be solved by a shift of any energy value up to 
this value. 

Substituting (37) in (25) and making use from: 



n z 

n=l 



Poly\og(2,z) = J2 

n=l 

finally we obtain next function for the density of grand thermodynamic potential: 

1 2 d 

u(p,R,0) = - 2 Polylog{2,z)Y,aAR){-Q-^) j Ci (40) 

p j=o m 

Now we calculate the density of the average of the energy .We can label it as pa which is: 

E r gjeje 



2 Note that V is disappeared since we calculate density. 
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After substituting (37) in it and doing the integration we obtain: 

C 

Pa = ~^Polylog(2,z) 
The mean value of the number of the particles in the unit volume of system is: 



(41) 




(42) 



Since < z < 1 the value of n is not negative. Using s = y = — 1£ we have: 



d 



(43) 



dm 2 



3=0 



Because the mass dependence of g(e) is undetermined , so we can not conclude strongly 
about the sign of entropy. It can be negative or positive [14]. 

Negative entropy is problematic if we accept that the entropy is in the association with 
the measure of the number of microstates in statistical mechanics. The intuition of statistical 
mechanics requires that the entropy of all physical components to be positive. Besides if 
we consider the universe as a thermodynamic system, the total entropy of the universe 
including DE and DM should satisfy the second law of thermodynamics. The GSL for 
phantom and non-phantom DE has been explored in [12]. It was found that the GSL can be 
protected in the universe with DE. The GSL of the universe with DE has been investigated 
in [13,25,26,27] as well. 3 It was found that the equation of state of the DE is uniquely 
determined and the phantom entropy is negative [13]. For the phantom, we either run into 
negative entropy problem or the GSL is violated. 4 Considering the apparent horizon as the 
physical boundary of the universe, it was found that both the temperature and entropy can 
be positive for the DE, including phantom. But in the general framework that we discussed, 
where the cosmological metric has a slow varying time dependency, we can assume that 
the time contiguity of volume and temperature is venial and for this reason we obtained 
a negative entropy. We must note here that if in an ensemble of particles the volume and 
temperature were non static, the statistical behavior of this system is not in the context of 
the equilibrium thermodynamics. 5 



We have derived a grand canonical partition function for a DE's like EOS. The thermo- 
dynamic equations are generally expressed in terms of three variables: the temperature T , 



3 We may think the DE temperature is equal or proportional to the horizon temperature Tjj- 

4 It is more realistic to consider that the physical volume and the temperature of the universe are related, 
since in the general situation they both depend on the scale factor a(t).It was found that the apparent 
horizon is a good boundary for keeping thermodynamic laws [13]. 

5 As in a FRW standard model of cosmology, since our background metric is not static and thereupon 
the apparent horizon is time variable, then construction of a local quantum statistical and counting all 
microstates must affect the negative entropy and we can deduce a statistical proof for non negativity and 
non violation of GSL like Gong and collaborations proof [6] . 



VI: CONCLUSIONS 
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chemical potential \i associating with a conserved particle number N, and scalar curvature 
R. We also reviewed some known results for an ideal Bose gas in curved spacetime. Com- 
paring our result for a negative pressure matter with an ideal gas we concluded that there 
is a significance generic difference between the statistical distribution of negative pressure 
matters and ideal gases. Specifically , the dispersion relation for a dark energy like EOS is 
non quadratic (in contradiction with a non interactive ideal Bose gas) and comparing this 
dispersion function with the statistical model of SQM we observe that also in SQM systems 
this dispersion relation is valid in some physically accessible region of the system[19]. 

Also we have derived the mean value of the energy of the system in terms of the curvature 
. We show that if a matter field (Boson) obeys an EOS like Dark energy and in cosmological 
language violates the well known energy conditions in General relativity , the degeneracy 
function of this gas must be an implicit constant function of energy This function implicitly 
is dependent on the mass of the particles which is not a thermodynamic quantity and remains 
as an undetermined function in this formula. 

The derived thermodynamic equations has show that the entropy may be positive or 
negative. As it was stated in [6], this negativity of the entropy is hidden behind this signif- 
icance assumption that volume and temperature are independent. But we know that in an 
expanding FRW universe both are functions of time. Hence if we assume that our metric 
is quasi static or has a negligible dynamics, any phantom field of DE with an intelligible 
statistical testimony has a negative entropy. If we want to generalize our argument to a 
realistic cosmological case, we must work in a non equilibrium statistical regime. 
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